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Experiments  in  Objective  Aviation  Weather 
Forecasting  Using  Upper-Level  Steering 


l.  INTRODUCTION 

When  winds-aloft  data  became  routinely  available  to  meteorologists,  the  fore¬ 
casters  quickly  noted  the  similarity  between  the  motion  of  storm  systems  and  the 
air  flow  at  levels  of  3-6  km  (500-700  mb).  This  relationship  became  the  basis  for 
12-  to  24-hour  forecasting  using  "upper-level  steering."  This  practice  received 
physical  justification  when  barotropic  theory  was  developed,  particularly  when 
Fjortoftl  applied  the  theory  to  forecasting  in  the  form  of  a  graphical  500-mb  fore¬ 
cast  technique,  using  a  space-averaged  flow  for  steering  vorticity  parcels  (or  "ad¬ 
verting"  the  parcels).  Soon  after,  Estoque2  and  Reed3  developed  graphical  1000- 
mb  forecast  techniques,  applying  baroclinic  theory,  and  these  techniques  steered 
thermal  patterns  using  upper-level  wind  flow. 

During  the  1970s,  the  fine-mesh  numerical  forecast  products  improved  so  that 
they  now  serve  as  the  basis  for  virtually  all  prognostic  charts,  and  the  steering 


(Received  for  publication  7  December  1983) 

1.  Fjortoft,  R.  (1952)  On  a  numerical  method  of  integrating  the  barotropic  vorti¬ 

city  equation,  Tellus  4:179-194. 

2.  Estoque,  M.A.  (1957)  A  graphical  integration  of  a  two-level  model,  J.  Atmos. 

Sci.  14:38-42. 

3.  Reed,  R.  J.  (1980)  On  the  practical  use  of  graphical  prediction  methods,  Mon. 

Wea.  Rev.  88:209-218. 
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techniques  have  become  obsolete.  Unfortunately,  even  the  latest  fine-mesh  numei — 
ical  model  uses  a  grid  of  about  190  km.  Consequently,  t He  weather  disturbances  in 
the  model  are  fairly  large,  producing  gradual  changes  over  periods  of  15-30  hours 
or  more.  This  has  left  the  local  forecaster  on  his  own  to  predict  the  more  rapid 
weather  changes  that  occur  over  periods  of  several  hours.  Many  of  these  rapid 
changes  are  associated  with  the  passage  of  surface  fronts  and  troughs  or  with  the 
arrival  of  clouds  and  precipitation,  phenomena  that,  in  mid-latitudes,  are  closely 
related  to  the  traveling  cyclones  and  anticyclones.  Since  the  upper-level  steering 
principle  once  served  well  in  forecasting  the  motion  of  cyclones  and  anticyclones, 
a  natural  deduction  would  be  that  steering  would  also  work  for  the  finer  scale  "sub¬ 
synoptic"  scale  disturbances  associated  with  the  cyclones  and  anticyclones. 

In  fact,  some  indication  exists  that  upper-level  steering  can  be  used  with  small 
sr.ale  disturbances.  In  a  forecast  test,  Muench^  used  objectively  determined  mo¬ 
tion  vectors  as  well  as  700-mb  and  500-mb  winds  to  make  objective  short-range 
forecasts  of  cloudiness  and  precipitation,  using  geosynchronous  satellite  imagery. 
The  test  showed  marked  improvement  over  persistence  for  forecasts  more  than  2 
hours,  with  nearly  equal  skill  for  the  objective  motion  vector  and  the  two  upper- 
level  winds.  When  the  size  of  the  cloud  and  precipitation  patterns  in  typical  satel¬ 
lite  images  and  radar  displays  is  taken  into  account,  the  test  must  be  considered 
as  at  least  partial  verification  of  the  concept  of  upper-level  steering  for  small-  or 
"meso-"scale  disturbances. 

While  the  satellite  data  are  nearly  ideal  for  cloud  cover  prediction  and  the  IR- 
visible  channel  combination  can  specify  precipitation  areas  fairly  well,  a  number 
of  other  weather  parameters  that  cannot  be  easily  related  to  satellite  measure¬ 
ments  are  of  interest  to  aviation  forecasters.  In  particular,  ceiling,  visibility, 
and  surface  wind  are  very  important  to  operations,  and  forecasts  of  altimeter  set¬ 
ting,  temperature,  and  humidity  are  sometimes  needed.  To  use  steering  for  ob¬ 
jective  forecasts  of  these  parameters,  one  would  have  to  use  the  hourly  aviation 
weather  reports  to  define  the  patterns  and  the  radiosonde  data  to  define  the  upper- 
level  winds. 

To  test  the  feasibility  of  using  the  steering  concept  for  objective  short-range 
terminal  forecasts,  several  forecast  models  of  increasing  complexity  were  devel¬ 
oped  and  tested.  This  report  will  describe  these  steering  or  advection  models,  the 
tests  that  were  performed,  and  the  conclusions  that  can  be  drawn  from  the  tests. 


4.  Muench,  H.S.  (1981)  Short-Range  Forecasting  of  Cloudiness  and  Precipitation 
Through  Extrapolation  of  GOES  Imagery.  AFGL-TR-81-0218,  ~ 


2.  ADVECTION  EXPERIMENT  I 


By  the  fall  of  1981,  the  AFGL  McIDAS  facility  had  the  means  to  make  tests  of 
the  steering  or  advection  technique  for  surface  weather  parameters.  Hourly  wea¬ 
ther  observations  and  radiosonde  ascents  have  been  routinely  saved  on  disk  stor¬ 
age  and  dumped  to  magnetic  tape  twice  a  week.  Software  developed  by  the  Univer¬ 
sity  of  Wisconsin  was  installed  that  performed  objective  analyses  of  almost  any 
parameter,  developing  grid-point  values  from  irregularly  spaced  observations.  In 
addition,  a  routine  had  been  written^  that  computed  locations  of  upstream  trajec¬ 
tories  of  upper-level  winds,  at  1-hour  intervals,  and  that  sought  out  the  corre¬ 
sponding  weather  parameter  from  another  grid  field  to  produce  advection  forecasts. 

In  mathematical  terms,  if  a  quantity  X  is  conservative  for  some  given  hori¬ 
zontal  flow  field,  we  can  write 


dX 
dt  * 


0 


which  becomes 


ax ,  ax 
aT*  -v*  aS- 


(i) 


where  V*  is  the  advecting  wind  speed  along  the  streamline  S  (or  steering  current). 
If  Eq.  (1)  is  integrated  from  time  t  =  0  to  t  =  t,  starting  at  location  x  =  0,  y  =  0.  we 

have  XQOt  -  X00Q  =  Xxyt  -  XQ00  or  XQQt  =  \^0.  which  means  that  the  forecast  for 
x  =  0.  y  =  0.  t  »  t  is  simply  the  upstream  value  of  X  at  x  =  x,  y  =  y,  t  =  0  (or  the 
initial  time). 

Using  the  same  notation,  a  forecast  of  "persistence"  means 

=  0  or  XQOt  »  XqoO  (2) 


A  3-week  test  of  the  advection  forecast  technique  was  conducted  from  mid -Oc¬ 
tober  to  early  November  1981,  using  the  McIDAS  facility.  Three  stations  were 
chosen:  Boston,  Mass.;  Burlington,  Vt.;  and  Atlantic  City,  N.  J.  Forecasts  were 
made  from  00  GMT  radiosonde  and  03  GMT  surface  data,  and  again  from  12  GMT 
radiosonde  and  15  GMT  surface  data.  Advection  forecasts  were  made  for  the  u  and 


5.  Wash,  C.H. ,  and  Whittaker,  T.M.  (1980)  Subsynoptic  analysis  and  forecasting 
with  an  interactive  computer  system.  Bull.  Am.  Meteorol.  Soc. 

81:1584-1591. 


v  components  of  the  surface  wind,  the  total  cloud  cover  (including  "thin"),  visibi¬ 
lity,  and  dewpoint.  While  the  McIDAS  routine  made  hourly  forecasts  out  to  15 
hours  in  1-hour  steps,  only  the  0-,  1-,  2-,  3-,  9-,  and  15-hour  forecasts  were 
recorded,  providing  adequate  time  coverage  and  simple  verification. 

a 

In  addition  to  persistence,  model-output-statistics  (MOS)  forecasts0  were  also 
used  as  a  control.  The  corresponding  MOS  forecasts  were  collected  from  the  Air 
Force  COMERS  circuit  for  the  3-,  9-,  and  15-hour  periods.  These  MOS  forecasts 
have  compared  favorably  with  local  subjective  forecasts  of  12-36  hours.  '  The 
forecasts  are  based  on  data  from  the  NMC  I.FM  model  prediction  together  with 
data  from  a  local  observation  taken  3  hours  after  the  initial  data  for  the  I.FM.  Sta¬ 
tistically  derived  relationships  combine  the  data  into  forecasts  of  a  wide  variety  of 
weather  parameters. 

In  the  verifications  for  this  study,  the  standard  of  comparison  is  the  root- 
mean-square  error  (rms  error).  The  rms  error  was  chosen  rather  than  the  mean- 
absolute-error  (MAE)  on  the  presumption  that  emphasis  should  be  given  to  reduc¬ 
ing  the  largest  errors.  In  general,  the  rms  errors  are  about  25-30  percent  larger 
than  MAEs.  In  some  cases  when  the  skill  scores  are  within  ±10  percent  of  zero, 
the  scores  based  on  rms  error  may  be  of  opposite  sign  to  those  based  on  MAE,  but 
normally  the  skill  scores  are  similar.  The  cloud  amount  and  visibility  scores 
were  based  on  the  MOS  categories  and  thus  represent  rms  category  errors.  The 
scores  for  wind  were  the  rms  vector  error  (combining  u  and  v  component  errors), 
and  the  dewpoint  scores  were  simply  rms  errors.  More  details  of  this  experiment 
may  be  found  in  the  Appendix  of  Muench.^  Two  different  advecting  wind  fields  were 
used,  the  700-mb  wind  and  the  850-300-mb  vertically  averaged  wind. 

Results  of  Advection  Experiment  I  are  summarized  in  Table  1,  and  show  the 
rms  errors  for  persistence  and  the  skill  scores  relative  to  persistence  [skill  = 
(persistence  error-forecast  error)/persistence  error].  Overall,  the  results  were 
somewhat  disappointing  and  difficult  to  understand.  For  wind,  cloud  amount,  and 
dewpoint,  the  only  positive  scores  (underlined  values),  representing  improvement 
over  persistence,  were  for  periods  beyond  3  hours,  and  even  these  were  less  ac¬ 
curate  than  the  MOS  forecasts.  If  the  steering  field  were  changing  with  time  or  if 
the  weather  patterns  were  undergoing  significant  development  or  decay,  then  the 
shorter  period  forecasts  should  have  been  more  accurate  than  the  longer  period 
forecasts;  instead,  we  see  the  reverse  here.  Only  the  visibilitv  forecasts  for  2  and 

6.  Glahn,  H.R.,  and  Lowry,  D.  (1972)  The  use  of  model  output  statistics  (MOS) 

in  objective  weather  forecasting,  J.  Appl.  Meteorol.  11:1203-1211. 

7.  Zurndorfer,  E.  A. ;  Bocchieri,  J.R.;  Carter,  G.  M.  ;  Dallaville,  J.  P.  ;  Gil  - 

hausen,  D.  B. ;  Hebenstreit,  K.  F. ;  Vercelli,  D.J.  (1979)  Trends  in  compar¬ 
ative  verification  scores  for  guidance  and  local  aviation  weather  forecasts, 
Mon.  Wea.  Rev.  107:799-811. 


Table  1.  Advection  Forecast  Experiment  I,  Skill  Scores  Relative  to  Persistence 


Wind  (rms  vector  error) 

Technique /Time  (h) 

0 

1 

2 

3 

9 

15 

Pereletence 

(10.  00) 

(±1.88) 

(±2.72) 

(±3.  18) 

(±3.  89) 

(±8.  11) 

Advection  700  mb 

(£2. 39 

-0.60 

-0.18 

-0.08 

+0.  15 

+0.21 

Advection  850-300  mb 

(±2.  39) 

-0.86 

-0.20 

-0.  09 

+0.  12 

+0.20 

AWS-MOS 

XXX 

XXX 

XXX 

+0.05 

+0.31 

+0.29 

Cloud  Cover  (rms  category  error,  0-3) 

Pereletence 

(to. 00) 

(±0.  38) 

(±0.  50) 

(±0.81) 

(±1.  ID 

(±1.22) 

Advection  700  mb 

(*0.73) 

-1.17 

-0.  84 

-0.  84 

+0.08 

-0.08 

Advection  850-300 mb 

(±0.73) 

-1.14 

-0.  80 

-0.  69 

+0.  00 

-0. 16 

AWS-MOS 

XXX 

XXX 

XXX 

-0. 19 

+0.  J 

>•0.27 

Visibility  (rms  category  error. 

0-4) 

Pereletence 

(±0.  00) 

f±0. 47) 

(±0.76) 

(±0.83) 

(±1. 

(±1. 18) 

Advection  700  mb 

(±0.48) 

-0.45 

+0.08 

+0.00 

-0. 

+0.  06 

Advection  850-300 mb 

(±0.48) 

-0.38 

+0.01 

+0.02 

-0.  1 

+0.09 

AWS-MOS 

XXX 

XXX 

XXX 

+0. 10 

+0.00 

+0. 11 

Dewpoint  Temperature  (rms  error. 

Celsius) 

Pereletence 

(±0.00) 

(to.  70) 

(±1.55) 

(±2.18) 

(±3.45) 

(±4.  93) 

Advection  700  mb 

(±1.61) 

-1.74 

-0.41 

-0.16 

-0.  08 

+0.  03 

Advection  850-300 mb 

(±1.61) 

-1.80 

-0.48 

-0.23 

-0.31 

+0.06 

AWS-MOS 

xxx 

XXX 

XXX 

XXX 

XXX 

XXX 

3  hours  show  promise,  and  these  skill  levels  are  less  than  MOS.  The  positive  skill 
scores  at  the  longer  times  indicate  that  some  features  are  being  properly  advected, 
and  the  negative  scores  at  shorter  periods  must  mean  that  some  features  of  small 
size  (or  short  lifetime)  are  being  improperly  advected,  wiping  out  some  gains  (al¬ 
beit  small)  over  persistence. 

There  are  three  principal  ways  that  anomalous  features  might  appear  in  the 
analyzed  fields,  leading  to  erroneous  forecasts.  First,  there  are  the  inevitable 
observational  errors  in  making  judgments  about  visibility  and  cloud  cover,  or  in 
determining  winds  and  dewpoints  from  strip-charts  or  dials.  These  errors  are 
likely  to  be  as  high  as  60  percent  of  the  rms  1-hour  change  or  persistence  error. 

Another  way  is  smoothing  some  features.  The  analysis  procedure  consists  of 
a  single  Cressman-type  grid-point  analysis®  on  a  1 -degree  longitude -latitude  grid 

8.  Cressman,  G.  P.  (1959)  An  operational  objective  analysis  system,  Mon,  Wea. 

Rev.  87:367-371. 
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(usually  15  by  33  degrees)  and,  admittedly,  does  oversmooth  some  important 
small-scale  features  'see  Gerlach^).  This  smoothing  accounts  for  many  of  the 
zero  hour  rms  errors  seen  in  parentheses  in  Table  1.  These  errors  represent 
the  inability  of  the  interpolation  routine  to  recover  the  initial  values  of  the  obser¬ 
vations  from  the  grid-point  values.  While  one  would  certainly  prefer  better  anal¬ 
yses,  perhaps  a  three-pass  Barnes-type  analysis^®  on  1/2-degree  grids,  such  a 
program  would  have  overtaxed  the  computational  capabilities  of  the  McIDAS  if  per¬ 
formed  on  the  15  by  33-degree  area  necessary  for  trajectories  out  to  15  hours. 

The  third  potential  source  of  error  is  the  presence  of  standing  or  stationary 
disturbances  in  the  otherwise  moving  weather  patterns.  Orographically  related 
cloud  patterns  were  known  to  have  degraded  some  of  the  forecasts  in  the  satellite - 
based  forecast  of  Muench,  ^  and  forecasters  are  quite  familiar  with  such  features 
as  lee-side  wind  (and  pressure)  troughs,  coastal  convergence  zones,  and  system¬ 
atic  lowering  of  visibility  with  altitude  (in  poor  weather).  Any  forecast  procedure 
that  attempted  to  advect  these  stationary  disturbances  would  certainly  suffer  from 
increased  forecast  errors. 


3.  ADVECTION  EXPERIMENT  1 1 

The  first  modification  of  the  advection  routine  was  an  attempt  to  reduce  errors 
incurred  when  stationary  disturbances  were  advected.  To  begin,  the  predictand  X 
is  now  assumed  to  be  composed  of  two  components,  a  stationary  or  orographic 

/N  /\ 

component  X  and  a  moving  component  X'.  Thus  X  =  X  +  X'.  The  local  time  deri¬ 
vative  is 


x  §2L.  ax: 

t  '  at  at 


We  now  assume  that  ax/at  =  0  (or  at  least  is  negligible  compared  with  dX'/dt). 
As  before,  we  assume  the  moving  pattern  (now  X’)  is  conservative  or  dX'/dt  =  0. 
This  also  means  that  d/dt  <aX'/a t)  =  0.  Since  SX/S t  =  0  everywhere  and  for  all 
time,  we  can  write  d/dt  ( dX/dt )  =  0.  Thus 


/ax\  a  (dx>  ax^ 
(eTf  dt  \  WT  +  ~ ; 


9.  Gerlach.  A.,  ed.  (1982)  Objective  Analysis  and  Prediction  Techniques. 
AFGL-TR-82-0394,  A"B"A  l3  - - 

10.  Barnes,  S.  L.  (1964)  A  technique  for  maximizing  details  in  numerical  weather 
map  analyses,  J.  Appl.  Meteor.  3:396-409. 


This  now  means  that  (dX/dX)  is  conservative,  and  we  can  write 


JLtex'v  „„  d_/ax\ 
d\  \at  )  w  ds\W) 

Integrating  over  time  (advecting  dX/di),  we  have 


A  second  integration  and  converting  an  integration  to  a  finite  difference  summation 
yields 

^00t  ^000  ZT  S^xyO  ^xy-Atj 

The  new  forecast  equation  is  then 

X00t  ■  X000  +irS(X,jrO  •  XJcy-it)  (4» 

that  we  will  call  the  "change-advection"  forecast  technique.  Following  the  lead  of 
early  synoptic  meteorologists,  who  adopted  3  hours  to  compute  pressure  tendency, 
we  chose  3  hours  for  At.  This  forecast  equation  has  the  desirable  feature  of  start¬ 
ing  with  the  observed  weather  at  the  forecast  site  at  time  zero,  thus  incorporating 
some  local  effects  and  starting  on  an  equal  footing  with  persistence. 

A  feasibility  test  of  the  change-advection  technique  was  conducted  for  12  cases 
from  early  December  1982  through  March  1983.  The  cases  were  not  random; 
nearly  all  involved  major  storm  systems  in  the  eastern  United  States  that  had  ap¬ 
preciable  motion  and  that  were  likely  to  contain  mixtures  of  stationary  and  moving 
patterns.  In  this  test,  forecasts  were  again  made  for  surface  wind  and  cloud  cover. 
However,  this  time,  temperature  replaced  dewpoint  for  a  test  of  how  both  advection 
and  change-advection  scored  on  a  parameter  with  an  important  diurnal  component. 
Again,  700-mb  flow  was  used  for  an  advective  wind,  but  now  the  700-500-mb  ver¬ 
tically  averaged  wind  was  given  a  spatial  smoothing  through  several  passes  of  a 
simple  5-point  averaging  routine  (called  "filter"  on  McIDAS).  In  the  previous  test, 
frontal  zones  had  not  been  forecast  eastward  fast  enough,  probably  because  of 
changing  steering  currents,  and  it  was  hoped  that  this  smoothing  (analogous  to  the 
Fjortoft  2  field)  might  correct  the  problem.  It  was  felt  that  700-mb  and  500-mb 
winds  provided  sufficient  information  on  vertical  structure  and  would  produce  a 
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vertical  average  wind  nearly  identical  to  the  850-300-mb  mean  used  in  Experi¬ 
ment  I. 

For  this  test,  the  number  of  stations  was  increased  to  six,  including  Boston, 
Mass. ;  Philadelphia,  Pa, ;  Washington,  D.  C. ,  National  Airport;  Buffalo,  N.Y.; 
Pittsburgh,  Pa.,  and  Columbus,  Ohio.  Four  advection  forecasts  were  made: 

700  mb  advection,  700-500  mb  (filtered)  advection,  700  mb  change -advection,  and 
700-500  mb  (filtered)  change-advection.  Persistence  and  the  NMC  MOS  forecasts 
(then  available  on  McIDAS  through  new  software®)  were  controls.  Forecasts  were 
made  for  0,  1,  2,  3,  6,  9,  12,  and  15  hours,  with  the  6-  and  12-hour  forecasts 
added  to  better  define  where  the  positive  skill  scores  began.  Six  cases  were  based 
on  00-03  GMT  data,  and  six  on  12-15  GMT  data.  Again,  rms  errors  were  com¬ 
puted  to  use  as  a  basis  for  comparing  scores.  The  results  of  Experiment  II  are 
summarized  in  Table  2  (in  the  same  format  as  Table  1). 

Overall,  the  results  of  this  experiment  were  encouraging.  For  the  winds,  we 
note  that  there  were  modest  positive  skill-scores  for  simple  advection  starting  at 


Table  2,  Advection  Forecast  Experiment  II,  Skill  Scores  Relative  to  Persistence 


VECTOR  WIND  (rms  error,  mps.  In 

parentheses) 

TECHNIQUE/TIME  <H) 

0 

1 

2 

3 

6 

9 

12 

15 

Persistence 

(±0.00) 

(12.01) 

(12.43 

(12.71) 

(14.30) 

(15.21) 

(16.04) 

(±6.25) 

Advection  TOO  mb 

(±1.51) 

-0.01 

-0.  11 

10.  05 

10.23 

♦0. 32 

♦0.28 

♦0.31 

Advection  7-5F 

(11.51) 

-0. 10 

-0.09 

iQ.  03 

10.26 

♦0.31 

10.29 

40.  29 

Ch ng  Adv.  700  mb 

(±0. 00) 

10.00 

iQ.  00 

-0.  15 

-0.06 

40.  00 

10.  05 

40.  00 

Chng  Adv.  7-5F 

(±0. 00) 

+0.  00 

♦0. 03 

-0.  13 

-0.08 

-0.  07 

-0.  03 

-0.  11 

NMC-MOS 

(10. 00) 

XXX 

XXX 

-0.  05 

it  mr 

40.  44 

XXX 

40.  47 

TOTAL  CLOUD  AMOUNT  (rms  category  error. 

0-3) 

TECHNIQUE/TIME  (H) 

0 

1 

2 

3 

6 

9 

12 

15 

Persistence 

(10.00) 

(10.  50) 

(10.60) 

(10.  70) 

(11.00) 

(11.  10) 

(11.30) 

(11.50) 

Advection  700  mb 

(10.  55) 

-0.  60 

-0.  17 

-0.05 

-0. 15 

-0.  09 

-0.05 

+0.03 

Advection  7-5F 

(10.  55) 

-0.50 

-0.  17 

-0.  14 

-0. 10 

-0.  09 

40.  08 

40.00 

Chng  Adv.  700  mb 

(10. 00) 

10.  10 

10.08 

40.29 

10.  10 

-0.  05 

40.  00 

40.00 

Chng  Adv.  7-5F 

(10. 00) 

-0.  20 

10.  00 

40.  00 

♦  0.05 

40.00 

10.04 

-0.  03 

NMC-MOS 

(10. 00) 

XXX 

XXX 

-0.29 

XXX 

40.  09 

XXX 

40.  27 

TEMPERATURE  (rma  error 

Celsius, 

in  parentheses) 

TECHNIQUE/TIME  (H) 

0 

1 

2 

3 

6 

9 

12 

15 

Persistence 

(10.  00) 

(10.86) 

(11.45) 

(11.96) 

(±2.75) 

(±3. 17) 

(13.63) 

(±4.  33) 

Advection  700  mb 

(±0.  92) 

-0.  74 

‘0.  55 

-0.  51 

-0.  13 

-0.01 

-0.23 

-0.  32 

Advection  7-5F 

(ID.  92) 

-0.  74 

-0.  53 

-0.  48 

-0.35 

-0.  12 

-0.  15 

-0.  09 

Chng  Adv.  700  mb 

(lo.  00) 

40.  19 

+  0.  17 

10.  18 

-0.01 

-1.09 

-1.  10 

-1.26 

Chng  Adv.  7-5F 

(10.  00) 

10.  19 

10.  17 

10.20 

40.  04 

-0.48 

-1.08 

-1.25 

NMC-MOS 

(10.00) 

XXX 

XXX 

40.  15 

♦  0.20 

♦  0.32 

10.41 

♦  0.  48 
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3  hours,  and,  compared  to  Experiment  I,  the  scores  are  higher.  The  change-ad- 
vection  scores  are  better  than  simple  advection  at  1  and  2  hours,  but  worse  than 
advection  for  longer  periods.  One  should  note  that  the  forecast  skill  scores  for 
MOS  are  also  higher  for  this  test  than  for  Experiment  I,  suggesting  that  the  wind 
patterns  in  these  selected  cases  were  moving  more  predictably. 

The  cloud  amount  forecasts  for  advection  had  positive  skill  scores  only  for  the 
longest  time  periods,  consistent  with  the  Experiment  I  results.  However,  the 
change-advection  technique  shows  positive  skill  for  the  1 -4-hour  period,  at  least 
for  700  mb  steering.  Since  this  test  was  based  on  selected  cases  similar  to  those 
in  the  test  using  satellite  imagery,  ^  we  might  have  expected  to  see  advection  better 
than  persistence  beyond  2  hours.  The  relatively  poor  performance  of  simple  ad¬ 
vection  using  only  airways  observations  for  periods  less  than  12  hours  suggests 
there  is  a  great  advantage  in  using  satellite  imagery  to  forecast  cloud  amount  ob¬ 
jectively,  no  doubt  due  to  the  much  higher  spatial  resolution. 

The  skill  scores  for  temperature  indicate  features  that  might  well  have  been 
anticipated.  The  advection  forecasts  are  worse  than  persistence  for  all  time  per¬ 
iods,  with  much  of  the  error  resulting  from  attempts  to  advect  the  orographically 
reduced  temperatures  of  mountain  stations  within  the  Appalachian  region.  Also, 
some  of  the  stations  in  this  test  are  "heat  island"  stations,  with  local  anomalies 
that  do  not  move.  As  might  be  expected,  the  change  advection  skill  scores  are 
much  better  (at  least  for  1-3  hours)  and  certainly  look  encouraging.  Unfortunately, 
the  change-advection  forecasts  have  very  poor  skill  scores  beyond  6  hours  and  ap¬ 
pear  to  "blow  up. " 

An  immediate  suspicion  is  that  all  the  change-advection  routine  is  really  doing  . 
with  temperature  is  extrapolating  the  diurnal  change  of  the  previous  3  hours  (aver¬ 
aged  spatially)  out  to  15  hours,  hour  by  hour.  During  the  first  few  hours,  the  fore¬ 
casts  benefit  by  this  extrapolation  as  the  diurnal  trend  continues,  but,  before  long, 
the  real  trend  reverses  and  the  forecasts  are  useless. 

A  further  indication  of  the  effects  of  the  diurnal  temperature  change  can  be 
seen  in  Figures  la  and  lb.  These  figures  depict  the  rms  errors  and  the  mean- 
arithmetic-errors  (or  bias)  for  15-GMT  forecasts.  The  rapid  increase  in  change- 
advection  rms  error  beyond  6  hours  is  very  evident  in  Figure  la.  In  Figure  lb,  it 
is  equally  evident  that  nearly  all  of  this  increase  is  the  result  of  the  nearly  linear 
increase  in  the  bias.  In  fact,  the  relative  standings  of  all  three  forecast  techniques 
in  Figure  la  are  essentially  determined  by  the  relative  size  of  the  bias  errors  at 
all  time  periods.  One  should  note  that  the  diurnal  temperature  rise  (negative  per¬ 
sistence  bias)  from  15  GMT  to  18  GMT  was  only  +1.4°  C.  This  temperature  rise 
is  less  than  the  normal  for  this  region  in  winter  and  reflects  effects  of  the  exten¬ 
sive  cloud  cover  for  these  12  cases.  And,  of  course,  the  winter  diurnal  is  smaller 
than  the  diurnal  for  the  other  seasons,  when  solar  heating  is  stronger.  Neverthe- 


RMS  Forecast  Error  vs  Time  •  36  Forecasts  Starting  1500GMT 
Forecast  time  (hours) 
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Figure  la.  RMS  Temperature  Errors  vs  Forecast  Time,  Advection 
Forecast  Experiment  il 


less,  even  in  winter,  the  diurnal  heating  cycle  plays  a  very  important  role  in  the 
short-range  temperature  forecasts. 

A  last  point  is  that  the  skill  scores  in  Table  2  do  not,  overall,  indicate  any 
obvious  advantage  to  either  700  mb  steering  or  to  the  vertical -and-space-averaged 
flow  (7-5F).  The  differences  between  the  skill  scores  of  these  two  are  small  and 


Mean  -  Arithmetic  -  Error  •  BIAS  <C> 


Forecast  Bias  vs  Time  -  36  Forecasts  starting  isoocmt 
Forecast  Time  -  (hours) 


Figure  lb.  Temperature  Bias  vs  Forecast  Time,  Advection  Forecast 
Experiment  II 


do  not  seem  to  be  consistent  over  all  time  periods.  Based  on  these  results,  there 
is  little  to  choose  between  these  two. 


4.  ADVECTION  EXPERIMENT  1 1 1 

The  effects  of  the  diurnal  temperature  changes  on  the  accuracy  of  the  advec- 
tion  techniques  were  to  be  expected.  One  would  not  realistically  attempt  to  develop 
a  temperature  forecast  technique  without  allowing  for  the  diurnal.  Radiative  heat¬ 
ing  and  cooling  directly  cause  the  temperature  change;  the  temperature  change,  in 
turn,  affects  the  relative  humidity;  and  the  humidity,  in  turn,  affects  the  visibility 
and  the  cloud  ceiling  height.  The  temperature  change  also  affects  the  hydrostatic 
stability  that,  in  turn,  affects  turbulent  exchange  and  the  surface  wind  speed.  Thus, 
several  other  parameters  also  have  diurnal  changes,  even  though  the  changes  might 
not  be  as  marked  as  the  daily  temperature  cycle.  In  the  change-advection  forecasts, 
small  diurnal  changes  become  amplified  as  each  iteration  is  carried  out.  For 
these  reasons,  one  of  the  first  steps  in  developing  Advection  Model  III  was  to  re- 
derive  the  forecast  equations  to  include  diurnal  changes. 

To  start,  we  now  say  X  =  X  +  X  +  X',  where  the  diurnal  component  is  indicated 
by  X,  To  obtain  a  persistence  forecast  with  diurnal  modification  (persistence  + 
diurnal),  we  assume,  as  before,  dX/dt  =  0  and  dX'/dt  =  0,  and  we  are  left  with 
dX/dt  =  dX/dt.  Integrating  this  last  expression  over  time,  the  forecast  expres¬ 
sion  becomes 

X00t  =  X000  "  X000  +  X00t  (5) 

Simply  put,  the  persistence  +  diurnal  forecast  consists  of  the  present  value  minus 
the  present  diurnal  component  plus  the  future  diurnal  component. 

To  obtain  the  expression  for  advection  forecast  modified  for  diurnal  effect 
(advection  +  diurnal),  we  start  with  X  =  X  +  X'  (assuming  X  =  0).  Since  X'  is  de¬ 
fined  as  the  traveling  component  and  is  assumed  to  be  conservative,  we  see  that 
dX/dt  =  dX/dt  or  d(X  -  X)/dt  =  0. 
d(X  -  X)/at  =  -V*d(X  -  X)/ d s. 
the  following  forecast  equation; 

x00t  =  Xxy0  "  Xxy0  +  x00t 

The  addition  of  the  diurnal  term  does  complicate  the  change-advection  forecast 
equation.  To  derive  a  change-advection  +  diurnal  formula,  we  start  with  a  conserv- 


Since  (X  -  X)  is  conservative,  we  have 
Integrating  this  expression  over  time,  we  have 


(6) 


ative  equation  d/dt  (dX/dt  -  dX/dt)  =  0,  from  which  we  derive 


After  the  necessary  integrations  over  time,  the  forecast  equation  is  found  to  be 

•^OOt  =  ^000  ”  ^000  +£t"  ^  ^xyO  “  ^xyO  “  ^xy-At  +  ^xy-Atj  +  ^00t  (7) 

These  new  forecast  equations  reached  the  state  of  complexity  where  a  com¬ 
plete  test  using  the  McIDAS  facility  would  have  required  extensive  software  devel¬ 
opment  and  modification.  Thus,  a  new  procedure  was  established:  analyzed  grid- 
fields  for  upper  level  winds  and  surface  weather  parameters  were  computed  on 
McIDAS  and  then  placed  on  magnetic  tape  for  further  processing  on  the  AFGL  C'DC 
computers.  In  addition,  the  history  file  of  hourly  surface  reports  was  copied  onto 
magnetic  tape,  to  be  used  later  in  the  forecasts  and  verification.  This  new  proce¬ 
dure  was  more  fully  automated  and  enabled  forecasts  to  be  made  for  more  stations. 
The  software  used  for  decoding  the  tapes  has  potential  use  for  future  routines  such 
as  error  checking,  data  filtering,  and  use  with  large  scale,  high  resolution  objec¬ 
tive  analysis  programs. 

In  the  discussion  of  the  results  of  Experiment  II  (Section  3),  there  are  sugges¬ 
tions  that  some  of  the  improvement  in  the  advection  scores  over  the  first  test  re¬ 
sulted  from  the  choice  of  cases.  The  advection  procedure  may  very  well  score 
better  relative  to  persistence  on  large,  steadily  moving  storm  systems.  The  plan 
for  Experiment  III  was  to  test  the  forecast  techniques  in  a  more  routine  fashion, 
taking  on  all  weather  situations  as  they  arose. 

The  month  of  March  1983  was  selected  for  the  test  because  the  data  sets  were 
more  complete  for  that  month  than  those  gathered  during  the  preceding  months. 

Six  new  forecast  stations  were  added  to  those  of  Experiment  II:  Windsor  Locks, 
Conn. ;  Allentown,  Pa.;  Flint,  Mich.;  Cleveland,  Ohio;  Covington,  Ky. ;  and  Rome 
AFB,  N.Y.  Forecasts  were  computed  and  verified  for  all  hours  from  0-15  hours, 
even  though  most  MOS  forecasts  were  only  made  for  3,  9,  and  15  hours  (tempera¬ 
ture  and  dewpoint  are  every  3  hours  in  NMC-MOS).  The  forecast  parameters  in¬ 
cluded  vector  wind,  wind  speed,  cloud  amount,  visibility,  temperature,  and  dew¬ 
point.  Visibility  was  converted  to  the  natural  log  of  visibility  (LnV)  to  place  more 
emphasis  on  the  numerically  small  changes  (and  errors)  that  occur  at  the  low  end 
of  the  visibility  range  and  that  are  so  operationally  important.  Actually,  the  cate¬ 
gory  ranges  in  the  MOS  forecasts  for  visibility  are  approximately  logarithmic  in 
their  distribution.  In  addition  to  the  vector  wind,  the  wind  speed  was  added  as  a 


separate  forecast  parameter,  obtained  numerically  in  the  advection  routines  from 
the  magnitude  of  forecast  wind  components. 

Having  the  upper-level-wind  grids  available  to  the  AFGU  CDC  6600  computer 
allowed  for  developing  more  elaborate  space-smoothing  routines.  A  new  space- 
smoothed  steering  field  was  created  by  starting  with  McIDAS-generated  2-degree 
latitude -longitude  grids  of  700-mb  and  500-mb  u  and  v  components  of  the  wind 
(United  States  east  of  105  W).  Then  vertically  integrated  grids  were  computed  us¬ 
ing  60  percent  of  the  700-mb  wind  components  and  40  percent  of  the  500-mb  wind 
components  (an  approximate  density  weighting).  The  grids  were  then  given  an 
8- by  6-point  spatial  smoothing.  Finally,  1 -degree  u-  and  v-component  grid  fields 
were  interpolated  for  use  in  the  advection  forecasts.  This  new  steering  field  was 
designated  the  700-500S  (where  S  indicated  space-smoothed),  and  became  the  third 
steering  field  in  addition  to  the  700  mb  and  the  7-5F  (or  McIDAS  filtered  steering 
fields).  With  three  steering  fields,  simple  advection  and  change-advection,  with 
and  without  diurnal  modification,  there  were  12  different  advection  forecast  tech¬ 
niques  to  test. 

The  derived  forecast  Eqs.  (5),  (6),  and  (7)  were  based  on  the  assumption  that 
the  diurnal  variations  of  the  weather  parameters  would  be  specified.  The  task  then 
is  to  produce  algorithms  that  would  specify  the  diurnal  component  of  the  variable 
under  differing  weather  conditions.  In  particular,  one  looks  for  a  specification  in 

the  form  of  X  »  f(t)  *  g  (Xj,  X2. - Xn)  where  Xj  represents  the  other  variables 

that  would  affect  the  amplitude  of  the  diurnal  curve  (that  is,  cloud  cover  affecting 
temperature  diurnal  range).  Few  climatological  summaries  contain  data  based  on 
time  of  day,  but  some  information  was  found  in  the  following  sources:  Air  Force 
Revised  Uniform  Summary  of  Station  Weather  Observations  (RUSSWO);^*  ^ 
National  Weather  Service  Decennial  Census  of  the  United  States  Climate:  Summary 
of  Hourly  Observations;^  and  the  Canadian  Hourly  Data  Summaries. 14  The  only 
hourly  mean  monthly  dry-bulb  and  wet -bulb  temperatures  were  found  in  the  Cana¬ 
dian  Summaries,  while  hourly  cloud  frequency  (four  categories)  values  were  taken 
from  the  Decennial  Census.  Otherwise,  data  presented  in  these  summaries  is  in 
the  form  of  3-hour  mean  frequencies  for  a  number  of  categories  for  variables  such 
as  ceiling,  visibility,  wind  direction,  and  wind  speed. 


11.  Revised  Uniform  Summary  of  Surface  Weather  Observations,  A -F.  U.  S.  Air 

Force,  AWs  (by  station)  (about  500  pages  each,  available  through  NTIS). 

12.  An  Aid  for  Using  the  Revised  Uniform  Summary  of  Surface  Weather  Observa- 

tions  (RUSSWO's)  (1983)  USAF  TAC/TN-83/001.  Scott  A FB.  IL  62225. 

13.  Climatography  of  the  United  States  No,  82-41.  Decennial  Census  of  United 

StatesClimate:  Summary  of  Hourly  Observations  1951-1960  (bv  station). 
Superintendent  of  bocuments,  TITS’.  GPO,  Washington,  D.  C. 

14.  Hourly  Data  Summyies  No.  1-84  (by  station).  Climatology  Division  Meteoro- 

logical  branch.  Department  of  Transport,  Canada. 
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Mean  Temperature  (Celsius)  March 


Some  hourly  temperature  values  for  Canadian  airports^  were  extracted  and 
are  shown  in  Figure  2.  The  curves  are  fairly  smooth.  Therefore,  3-hour  means 
would  describe  the  overall  features  but  would  underestimate  the  amplitude  because 
of  the  averaging.  Also,  there  would  be  difficulties  positioning  the  time  of  maxi¬ 
mum  and  minimum. 

Because  of  the  lack  of  resolution  in  the  3 -hour  data  and  the  difficulties  work¬ 
ing  with  frequencies  (rather  than  mean  values),  and  because  the  climatological 
values  are  functions  of  time  only,  more  detailed  diurnal  information  was  developed 
by  using  the  MclDAS  archived  hourly  observations.  A  program  extracted  hourly 
reports  from  30  stations  east  of  90  W  from  the  archive  tape  for  March  1983. 
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Figure  2.  Mean  Temperature  vs  Time  of  Day,  4  Canadian  Stations,  March 


The  mean  values  for  1-hour  changes  in  wind  speed,  natural  log  of  visibility 
(LnV),  and  temperature  were  derived  for  a  variety  of  initial  weather  conditions, 
including  cloud  cover  (by  amount  and  layer),  wind  speed,  and  dewpoint  spread.  To 
compute  a  diurnal  curve,  the  first  step  was  to  subtract  the  24-hour  mean  of  the 
1-hour  changes  from  individual  changes  in  order  to  remove  any  bias.  For  example, 
clear-sky  cases  may  contain  cooling  due  to  advection.  Next,  a  first -guess  diurnal 
curve  was  constructed  by  integrating  the  1-hour  changes  (bias  removed)  starting 
at  OOGMT.  Because  OOGMT  is  rarely  a  null  point,  a  phase  correction  was  made 
by  computing  the  24-hour  average  and  subtracting  from  individual  values.  Figures 
3  and  4  show  examples  of  diurnal  curves  for  three  different  cloud  conditions  for 
temperature  and  wind  speed.  These  curves  represent  the  mean  of  all  80  stations. 

The  systematic  differences  in  diurnal  amplitude  going  from  to  low-broken/- 
overcast  are  quite  noticeable  in  Figures  3  and  4.  In  the  case  of  temperature,  the 
shapes  and  phase  for  all  curves  are  similar,  indicating  that  a  fairly  simple  algo- 


Figure  3.  Mean  Temperature  Departure  vs  Time  of  Day,  3  Cloud  Conditions.  March 
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rithm  of  the  form  T  =  f(t)  g(cloud)  could  be  constructed.  One  problem  is  that  the 
weather  elements  that  affect  diurnal  cycles  contain  diurnal  components  themselves, 
and  it  is  difficult  to  know  where  to  begin.  The  simplest  procedure  seemed  to  be  to 
start  with  the  cloud  cover.  Cloud  cover  has  a  small  diurnal  amplitude  compared  to 
its  dynamic  changes  and,  consequently,  neglect  of  diurnal  feedback  would  not  be 
serious.  The  following  algorithms  were  developed,  based  on  Decennial  Census 


hourly  cloud  amount  frequencies: 


C  =  C  -  C  =  0.  3  *  (30  -  C)  *  f(t)  (C  >  10) 
C  =  C  -  C  =  0.  6  *  C  *  f(t)  (C<10) 


where  C  is  the  non-diurnal  cloud  amount,  which  ranges  from  0  for  clear  to  30  for 
overcast.  If  C  and  not  C  is  known,  these  equations  can  easily  be  solved  to  obtain 
C  and  C. 

The  March  1983  study  based  on  80  stations  indicated  that  the  strongest  effect 
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Figure  4.  Mean  Wind  Speed  Departure  vs  Time  of  Day,  3  Cloud  Conditions,  March 
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on  diurnal  wind  speed  cycle  was  the  cloud  cover,  defined  above.  With  the  wind 
speed  S  given  in  tenths  of  knots,  the  derived  algorithm  is 

S  -  S  -  S  =  (43.  -  C)  *  f(t) 

In  Figure  4,  the  diurnal  curve  for  wind  speed  with  clear  sky  is  very  similar  to 
that  of  temperature,  implying  a  common  f(t).  However,  for  large  cloud  amounts, 
there  appears  to  be  a  phase  shift  in  the  peak  wind  speed.  In  this  experiment,  this 
phase  shift  was  ignored  because  the  errors  that  would  be  incurred  would  be  rela¬ 
tively  small  and  the  programming  would  be  considerably  simpler. 

The  March  1983  data  were  used  to  develop  the  following  algorithm  for  temper¬ 
ature,  where  T  is  expressed  in  tenths  of  degrees  Celsius: 

T  =  T  -  T  =  (5000  -  116  *  C)  *  f(t)/(l  +  7.  3  *  /T) 

(Values  of  3  less  than  10  or  1  knot  are  made  equal  to  10. ) 

Next  the  values  for  the  diurnal  V  =  10  *  LnV  are  approximated  by 

V  =  V  -  V  =  (31.  -  .046  *  V)  *T/?5 

If  V  and  not  V  is  known,  the  equation  can  still  be  easily  solved.  This  formulation 
is  probably  an  oversimplification,  and,  since  there  were  few  very  low  visibilities 
in  the  data  base,  one  should  not  expect  good  performance  with  dense  fog. 

Difficulties  were  encountered  in  trying  to  find  obvious  diurnal  patterns  in  wind 
direction.  Averaged  over  all  stations,  there  was  some  tendency  for  a  counter¬ 
clockwise  shift  of  a  few  degrees  in  the  evening  and  a  clockwise  shift  in  mid-morn¬ 
ing.  However,  the  amplitudes  were  of  the  order  of  4  degrees  for  strong  winds  and 
10  degrees  for  weak  winds.  Such  small  values  would  have  little  effect  on  the  vec¬ 
tor  wind  changes  from  hour  to  hour.  Undoubtedly,  at  some  individual  stations, 
there  are  much  larger  systematic  wind  direction  diurnal  patterns  associated  with 
local  terrain,  but  the  advection  forecast  equations  would  be  difficult  to  solve  with¬ 
out  a  universal  expression.  Thus,  the  diurnal  amplitude  used  for  wind  direction  A 

'■V'  _ 

was  A  =  A  -  A  =  0. 

Similarly,  there  were  also  difficulties  finding  consistent  diurnal  dewpoint  tem¬ 
perature  patterns.  Again,  local  effects  may  be  important,  and,  in  addition,  a 
variety  of  sensor  problems  may  be  involved.  A  more  pronounced  dewpoint  diurnal 
cycle  would  probably  be  seen  in  the  summer,  but,  for  this  March  experiment,  we 
assumed  for  dewpoint  H  that  H  =  H  -  H  =  0. 

A  careful  look  at  the  shapes  of  the  diurnal  curves  reveals  obvious  problems  to 
fitting  simple  sine  or  cosine  functions.  Physically,  it  is  possible  that  the  combina- 


tion  of  rising  sun,  changing  boundary  layer  depth,  and  outgoing  radiation  could 
produce  an  approximate  sinusoidal  pattern  from  dawn  to  late  afternoon.  However, 
by  sunset,  the  process  changes  to  one  of  net  outgoing  radiation  that  decreases  in 
magnitude  as  the  temperature  falls,  more  like  an  exponential  decay  process.  For 
this  reason,  essentially  two  algorithms  were  constructed,  one  for  day  and  one  for 
night.  In  the  following  expressions,  t  represents  the  time  in  hours  past  0600  local 
sun  time. 

0615  to  1715  LST  f(t)  =  .  60  *  cos  [(t  -  9.  0)  .  06 

1715  to  0615  LST  fit)  =  .98  *  exp  !  -  .  2  *  (t  -  11.  5)]  -  .  564 

Besides  inclusion  of  the  diurnal  change  expressions  in  the  forecast  equations 
and  the  addition  of  the  700-500  mb  space-smoothed  steering  field  (700-500S),  one 
more  change  was  made  in  the  forecast  procedure.  One  problem  perceived  in  the 
forecasting  of  surface  wind  was  the  presence  of  high  frequency  fluctuations  in  the 
nominal  1 -minute  averaged  wind  transmitted  every  hour.  ^  These  fluctuations  are 
caused  by  boundary  layer  turbulence  and  would  probably  introduce  noise  into  the 
analyses  and  errors  into  the  forecasts.  Unfortunately,  the  data  in  the  hourly  re¬ 
ports  are  not  sufficient  to  properly  filter  out  the  turbulent  contribution.  An  attempt 
was  made  in  this  direction  by  time  averaging  the  u  and  v  grid  fields  with  those  of 
the  hour  before.  The  latest  and  prior  observations  should  also  have  been  averaged 
for  use  in  the  change-advection  routine,  but  this  step  was  neglected  in  the  forecast 
program. 

The  data  processing  consisted  of  generating  some  26  analysis  grids  on  the 
AFGL  McIDAS  for  each  of  the  12  cases  and  placing  the  grids  on  tape.  On  a  case- 
by-case  basis,  forecasts  were  generated  and  verified  on  the  AFGL  computer  (with 
the  scores  saved  on  disk  files).  After  all  cases  were  processed,  scores  for  the 
complete  experiment  were  computed.  Having  scores  for  individual  cases  allowed 
for  a  later  check  into  systematic  differences  for  different  overall  weather  condi¬ 
tions.  The  resulting  skill  scores  are  presented  in  Tables  3a,  3b,  3c,  3d,  3e,  and 
3f.  The  rms  forecast  errors  are  shown  graphically  in  Figures  5a,  5b,  5c,  5d,  5e, 
and  5f  for  persistence  and  the  700-500S  advection  and  change-advection  forecasts 
as  well  as  the  NMC-MOS  forecasts.  Figures  6a,  6b,  6c,  6d,  and  6e  illustrate  the 
effects  of  adding  the  diurnal  algorithms  to  the  forecast  routines  by  examining  the 
mean-arithmetic  error,  or  bias. 

Overall,  the  results  of  Experiment  III  are  quite  similar  to  Experiments  I  and 
II.  The  advection  forecasts  for  vector  wind  show  positive  skill  beginning  about  2- 


Muench,  H.S.  (1982)  An  Appraisal  of  the  Short-Range  Forecast  Problem 
Using  Power  Spectra,  AFGL-TR-82-0353,  ADA129315. 
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tion  of  rising  sun,  changing  boundary  layer  depth,  and  outgoing  radiation  could 
produce  an  approximate  sinusoidal  pattern  from  dawn  to  late  afternoon.  However, 
by  sunset,  the  process  changes  to  one  of  net  outgoing  radiation  that  decreases  in 
magnitude  as  the  temperature  falls,  more  like  an  exponential  decay  process.  For 
this  reason,  essentially  two  algorithms  were  constructed,  one  for  day  and  one  for 
night.  In  the  following  expressions,  t  represents  the  time  in  hours  past  0600  local 
sun  time. 

0615  to  1715  LST  f(t)  =  .  60  *  cos  [(t  -  9.  0)  .  06 

1715  to  0615  LST  f(t)  =  .98  *  exp  [  -  .2  *  (t  -  11.5)]  -  .564 

Besides  inclusion  of  the  diurnal  change  expressions  in  the  forecast  equations 
and  the  addition  of  the  700-500  mb  space-smoothed  steering  field  (700-500S),  one 
more  change  was  made  in  the  forecast  procedure.  One  problem  perceived  in  the 
forecasting  of  surface  wind  was  the  presence  of  high  frequency  fluctuations  in  the 
nominal  1 -minute  averaged  wind  transmitted  every  hour.  14  These  fluctuations  are 
caused  by  boundary  layer  turbulence  and  would  probably  introduce  noise  into  the 
analyses  and  errors  into  the  forecasts.  Unfortunately,  the  data  in  the  hourly  re¬ 
ports  are  not  sufficient  to  properly  filter  out  the  turbulent  contribution.  An  attempt 
was  made  in  this  direction  by  time  averaging  the  u  and  v  grid  fields  with  those  of 
the  hour  before.  The  latest  and  prior  observations  should  also  have  been  averaged 
for  use  in  the  change -ad vection  routine,  but  this  step  was  neglected  in  the  forecast 
program. 

The  data  processing  consisted  of  generating  some  26  analysis  grids  on  the 
AFGL  McIDAS  for  each  of  the  12  cases  and  placing  the  grids  on  tape.  On  a  case- 
by-case  basis,  forecasts  were  generated  and  verified  on  the  AFGL  computer  (with 
the  scores  saved  on  disk  files).  After  all  cases  were  processed,  scores  for  the 
complete  experiment  were  computed.  Having  scores  for  individual  cases  allowed 
for  a  later  check  into  systematic  differences  for  different  overall  weather  condi¬ 
tions.  The  resulting  skill  scores  are  presented  in  Tables  3a,  3b,  3c,  3d,  3e,  and 
3f.  The  rms  forecast  errors  are  shown  graphically  in  Figures  5a,  5b,  5c,  5d,  5e, 
and  5f  for  persistence  and  the  700-500S  advection  and  change-advection  forecasts 
as  well  as  the  NMC-MOS  forecasts.  Figures  6a,  6b,  6c,  6d,  and  6e  illustrate  the 
effects  of  adding  the  diurnal  algorithms  to  the  forecast  routines  by  examining  the 
mean-arithmetic  error,  or  bias. 

Overall,  the  results  of  Experiment  III  are  quite  similar  to  Experiments  I  and 
II.  The  advection  forecasts  for  vector  wind  show  positive  skill  beginning  about  2- 


15.  Muench,  H.S.  (1982)  An  Appraisal  of  the  Short-Range  Forecast  Problem 
Using  Power  Spectra,  AFGL-TR-82-0353.  ADA129315. 
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Figure  5a.  RMS  Wind  Vector  vs  Forecast  Time,  Advection  Forecast  Experiment  III,  700-500S  Flow 


Forecast  Interval  —  (  hours) 


Figure  5c.  R MS  Cloud  Amount  vs  Forecast  Time,  Advection  Forecast  Experiment  III,  700-500S  Flow 


RMS  L.n  Visibility  vs  Forecast  Time,  Advection  Forecast  Experiment  III,  700-500S  Flow 


Forecast  Interval  —  (  hour ) 


Figure  5e.  RMS  Temperature  vs  Forecast  Time,  Advection  Forecast  Experiment  III,  700-500S  Flow 


S  Dewpoint  vs  Forecast  Time,  Advection  Forecast  Experiment  III,  700-500S  FI 
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Figure  6b.  Persistence  Bias  vs  Forecast  Time,  Cloud  Amount,  Advection  Fore¬ 
cast  Experiment  III. 


biases  were  of  sign  opposite  to  that  of  simple  persistence  for  both  03 -GMT  and  15- 
GMT  forecasts.  The  overcompensation  appears  to  be  the  result  of  an  encoding  er¬ 
ror;  the  error  will  be  corrected  for  future  tests.  However,  considering  that  the 
diurnal  changes  in  speed  03  GMT  (or  15  GMT)  for  a  mixture  of  clear  and  cloudy 
cases  would  be  about  ±0.6  mps,  representing  a  small  part  of  persistence  error, 
even  a  correct  diurnal  forecast  probably  would  not  have  produced  more  than  a  few 
percent  increase  in  the  skill  scores.  (Of  course,  producing  even  a  few  percent 
increase  isn't  always  easy. ) 

In  Experiments  II  and  III,  wind  speed  itself  was  not  forecast,  so  the  relatively 
low  skill  scores  in  Table  3b  are  somewhat  surprising,  if  not  actually  puzzling. 

For  instance,  why  should  the  advection  scores  for  3  hours  be  positive  for  the  wind 
vector  but  negative  for  the  wind  speed?  The  inference  is  that  the  wind  direction 
changes  are  being  forecast  better  than  the  speed  changes  are. 

The  main  problem  appears  to  be  that  the  wind  speeds  observed  by  airways  sta¬ 
tions  are  not  as  compatible  with  each  other  as  the  wind  directions  are.  First,  not 


Bias  vs  Time  (hours) 


Figure  6c.  Persistence  Bias  vs  Forecast  Time,  Ln  Visibility,  Advection  Forecast 
Experiment  III 

all  sensors  are  mounted  at  the  same  height  above  ground.  This  leads  to  problems 
because  of  the  systematic  increase  in  speed  with  height.  A  second  difficulty,  de¬ 
scribed  in  some  detail  by  Fujita  and  Wakimoto, 16  is  that  the  effects  of  surface 
roughness  and  nearby  obstructions  on  wind  observations  can  vary  considerably 
from  site  to  site.  To  analyze  surface  wind  for  making  objective  forecasts  or  dy¬ 
namic  computations  such  as  divergence,  station -by -station  adjustment  of  wind 
speed  must  be  made,  and,  in  some  cases,  the  speed  adjustments  may  be  direction- 
dependent.  If  the  wind  speed  at  a  given  site  were  systematically  slow  or  fast,  at 

16.  Fujita,  T.  T. ,  and  Wakimoto,  R,M.  (1982)  Effects  of  miso- and  meso-scale 
obstructions  on  the  PAM  winds  obtained  during  project  NIMROD,  J,  of 
AppI.  Meteorol,  21:840-858. 
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Figure  6d.  Persistence  Bias  vs  Forecast  Time,  Temperature,  Advection  Fore¬ 
cast  Experiment  III 

least  the  sign  of  the  changes  (though  not  the  magnitude)  would  be  correct.  We  note 
in  Table  3b  that  the  change -advection  scores  for  wind  speed  are  noticeably  better 
than  the  scores  for  advection  in  the  1 -4-hour  period,  whereas  the  advection  and 
change -advection  scores  were  similar  for  the  vector  winds  (Table  3a). 

As  in  previous  experiments,  the  advection  skill  scores  for  total  cloud  amount 
are  low,  mostly  negative  (Table  3c).  The  diurnal  algorithm  could  probably  be  im¬ 
proved  a  little,  as  both  persistence  and  persistence  +  diurnal  had  similar  biases  as 
seen  in  Figure  6c.  However,  the  contributions  of  the  diurnal  cycle  to  overall 
changes  in  March  are  so  relatively  small,  the  effort  would  not  be  justified.  A  bet¬ 
ter  way  to  improve  the  cloud  forecasts  would  be  to  use  satellite  imagery,  with 
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Figure  6e.  Persistence  Bias  vs  Forecast  Time,  Dewpoint,  Advection  Forecast 
Experiment  III 

much  higher  spatial  resolution  (and  perhaps  better  consistency),  to  specify  the 
cloud  fields  in  the  advection  routines. 

The  results  for  the  visibility  forecasts  are  encouraging  (see  Table  3d).  The 
skill  scores  for  advection  are  positive  beyond  2  hours,  though  the  level  of  skill  is 
modest.  The  diurnal  algorithms  definitely  improve  the  persistence  forecasts  out 
to  9  hours,  as  indicated  by  the  skill  scores  in  Table  3d  and  the  smaller  biases  in 
Figure  6c.  During  the  12  forecast  cases  of  March  1983,  most  of  the  reduced  visi¬ 
bility  was  either  directly  or  indirectly  associated  with  precipitation  (mostly  rain), 
so  the  positive  skills  may  indicate  success  in  advecting  precipitation  areas,  which 
had  been  found  by  Muench^  in  the  satellite -based  advection  experiments.  Low  vi¬ 
sibility  in  precipitation  cases  usually  accompanies  low  cloud-ceiling  heights,  so 
there  is  reason  to  expect  that  the  advection  techniques  would  also  find  modest  suc¬ 
cess  forecasting  ceiling. 

As  expected,  the  effects  of  modifying  the  forecast  techniques  for  the  diurnal 
change  were  quite  pronounced  in  the  case  of  temperature.  The  skill  scores  for 
persistence  diurnal  are  themselves  high  enough  to  be  considered  useful.  In  addi¬ 
tion,  the  bias  errors  in  Figure  6d  indicate  that  further  improvement  might  result 
from  increasing  the  magnitudes  of  the  diurnal  temperature  change.  The  simple 
advection  scores  are  as  poor  as  those  in  Experiment  II,  but  the  addition  of  the  di¬ 
urnal  cycle  improved  the  scores  to  the  point  where  there  are  some  positive  scores 
at  4-6  hours  out  into  the  forecast.  The  change-advection  scores  are  again  quite 
good  out  to  4  hours,  but  the  addition  of  the  diurnal  cycle  only  improved  the  scores 
beyond  4  hours.  This  result  is  consistent  with  the  earlier  suggestion  that  all  the 
skill  for  the  change-advection  in  the  first  4  hours  comes  from  the  extrapolation  of 
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the  diurnal:  the  diurnal  modification  only  subtracted  out  the  diurnal  in  one  place 
and  added  it  back  in  another. 

The  dewpoint  forecast  skill  scores  are  as  low  as  in  Experiment  I,  and  there 
are  obviously  some  problems.  The  rms  changes  for  1-6  hours  are  somewhat  lower 
than  for  temperature,  so  one  problem  could  be  that  the  uncertainties  in  the  mea¬ 
surements  are  causing  changes  as  large  as  the  real  changes  (or  larger).  However, 
the  failure  of  advection  to  produce  skillful  forecasts  for  either  temperature  or 
dewpoint  at  12  hours,  when  changes  are  larger,  raises  a  question  of  whether  the 
wrong  steering  was  used.  Physically,  rather  than  using  TOO  mb  or  500  mb,  the 
850-  or  even  1000-mb  level  may  be  more  reasonable.  In  these  experiments,  the 
higher  levels  were  used  in  the  belief  that  they  would  better  portray  the  combined 
effects  of  advection  and  the  modifying  effects  that  clouds  and  precipitation  have  on 
temperature  and  dewpoint. 

5.  CONCLUSIONS  AND  DISCUSSION 

The  advection  techniques  tested  thus  far  have  not  shown  the  level  of  skill  rel¬ 
ative  to  persistence  that  one  would  like  to  see  for  an  operational  forecast  tool  in 
routine  use.  However,  for  the  2-  to  6-hour  period,  the  techniques  could  provide 
useful  guidance  for  wind  and  visibility,  especially  since  the  MOS  scores  are  low 
for  these  parameters. 

When  Experiment  III  was  conducted,  forecasts  were  also  verified  on  a  case- 
by-case  basis  to  see  if  performance  was  weather  dependent.  Asa  simple  test, 
standard  deviations  were  computed  about  spatial  means,  for  vector  wind  (uv), 
cloud  cover  (C),  temperature  (T),  and  dewpoint  (H),  The  standard  deviations  were 
correlated  to  3-hour  persistence  error  (or  change).  The  values  computed  for  the 
12  cases  were  rw  =  +.  60,  rc  =  +.49,  rT  =  -.  24,  and  rH  =  +.  56. 

For  wind,  cloud  cover,  and  dewpoint,  these  correlations  indicate  a  dynamic 
link  between  gradients  and  local  change  (perhaps  advection)  and  imply  that  a  fair 
portion  of  the  change  may  be  forecastable.  In  the  case  of  temperature,  there  was 
a  -0.66  correlation  of  rms  3-hour  change  to  area,  average  cloud  cover  again  show¬ 
ing  importance  of  radiation.  When  the  improvements  over  persistence  for  advec¬ 
tion  were  compared  to  the  spatial  variability,  the  correlations  were  ruv  =  +0.  30, 
rC  "  +0. 30,  r*j>  *  -0.47,  and  rjj  =  +0. 16.  These  small  positive  correlations  (ex¬ 
cept  temperature)  show  the  advection  routine  performs  better  when  gradients  are 
large  ("active-looking"  situations). 

Several  steps  are  still  to  be  taken  before  the  validity  of  the  advection  concept 
can  really  be  tested.  As  a  start,  the  basic  airways  reports  must  be  edited  before 
being  used  in  objective  forecast  procedures.  In  addition,  the  parameters  must  be 


"normalized"  to  reduce  local  effects  (particularly  important  for  wind  speed).  Tem¬ 
perature,  visibility,  and,  to  some  extent,  dewpoint  are  not  conservative  for  mo¬ 
tions  up  and  down  terrain,  so  some  transformations  should  be  made,  analogous  to 
"potential  temperature,  "  before  they  are  used  in  advection  routines.  Next,  the 
analysis  routine  used  was  probably  too  crude,  and  a  more  sophisticated  "Barnes - 
type"  scheme, 10  using  multiple  passes,  should  be  tried.  And,  of  course,  if  satel¬ 
lite  imagery  were  introduced  into  the  procedure,  it  would  certainly  improve  the 
cloud  forecasting. 

The  situation  in  short-range  forecasting  today  may  well  be  similar  to  the  sta¬ 
tus  of  synoptic -scale  forecasting  in  the  early  1950s,  when  prognostic  charts  were 
prepared  manually  without  computer  guidance.  At  that  time,  there  were  consider¬ 
able  differences  in  skill  from  one  forecaster  to  another,  reflecting,  to  some  extent, 
their  relative  abilities  at  separating  signal  from  noise  in  the  basic  data.  When  ob¬ 
jective  forecast  procedures  became  available,  they  equalled  the  performance  of  the 
poorer  forecasters  at  first,  and,  when  improved  initialization  procedures  were 
introduced  into  the  dynamic  models,  the  skill  reached  that  of  the  best  forecasters. 
(The  biggest  advantage  of  the  dynamic  model  is  kindling  a  hemisphere  at  many  lev¬ 
els  in  the  time  the  synoptician  forecasts  two  levels  over  the  United  States. )  In  the 
mesoscale,  for  the  1-  to  4-hour  period,  skill  scores  of  about  0. 1  to  0. 4  seem  to  be 
par  at  present4  for  routine  subjective  forecasts,  and,  as  in  the  1950s,  individual 
forecasters  differ  considerably  in  skill. 17  IT  the  analogy  is  applicable,  all  of  the 
suggested  improvements  just  given  might  only  boost  the  skill  scores  to  the  +0. 1  to 
+0. 4  range.  Unless  the  future  mesoscale  dynamic  models  can  diagnose  disturb¬ 
ances  the  forecasters  are  missing,  even  these  models  may  only  achieve  skills  rel¬ 
ative  to  persistence  of  about  +0. 4  to  +0. 5,  meaning  that  half  the  short -period 
changes  will  still  be  missed.  However,  regular  forecasts  of  half  the  changes  would 
still  be  much  better  than  none. 
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